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SUMMARY 


A method is developed that determines the placement of an active control surface for maximum effectiveness 
in suppressing flutter. No specific control law is required by this method which is based on the aerodynamic energy 
concept. It is argued that the spanwise placement of the active controls should coincide with the locations where 
maximum energy per unit span is fed into the system. The method enables one to determine the distribution, over 
the different surfaces of the aircraft, of the energy input into the system as a result of the unstable fluttering mode. 
The method is illustrated using three numerical examples. 

INTRODUCTION 

All active flutter suppression systems require sensors to detect the movement of the lifting surface and to activate 
a control surface according to a synthesized control law. Most of the work performed to date (Mahesh and others, 
1981; Mukhopadhyay and others, 1981; Newsom, 1979; Nissim and Abel, 1978; Nissim and others, 1976) relate to 
the development of control laws based on predetermined locations of sensors and control surfaces.These locations 
of sensors and control surfaces are determined either arbitrarily, or by means of a trial and error procedure (Nissim 
and others, 1976). 

The aerodynamic energy approach developed in Nissim (1971) indicates that the sensors should be on the same 
streamwise strip where the activated control surface is placed, much in accord with Wykes’ identically located ac- 
celerometer force (ILAF) concept (Wykes and Mori, 1965). The aerodynamic energy concept also indicates (Nissim, 
1977, p. 25) that the best chordwise location of a sensor activating a trailing edge control surface is around the 65- 
percent chord location. The best chordwise location for a sensor activating a leading edge surface is shown to lie 
upstream of the wing (20-percent chord upstream of the leading edge), or alternatively, two sensors located along 
the same chord should be used. 

Recent work relating to mass-unbalance effects of control surfaces indicates that the spanwise location of the 
sensors within the control surface strip is related to the spanwise location of the center of gravity of the control 
surface and to the centroid of the strip’s surface. The determination of the best spanwise placement of an active 
control surface for flutter suppression using a methodical approach rather than by trial and error is the primary topic 
of this paper. 

The authors describe a method that determines the aforementioned spanwise placement of the activated control 
surface without resorting to any specific control law. The method is based on the aerodynamic energy concept 
whereby the activated control surface is placed at the location where most of the energy is fed into the unstable 
structure. The following analysis will assume, for sake of simplicity, that the generalized aerodynamic forces are 
obtained using a doublet-lattice program. This is done for purposes of illustration only, and the method can be 
applied to any form of calculation of the generalized aerodynamic forces. 

NOMENCLATURE 

o t - area of box* 

b semichord length 

h»/ 3/4 deflection at 3/4 -chord point of box »' due to deflection of mode j 


I 


hih/*/b 


hih/A 

hij 1/A deflection of 1 /4 -chord point of box » due to deflection of mode j 

k reduced frequency (usb/V) 

£ number of strips on the lifting surfaces 

m number of boxes 

n number of modes 

Q dynamic pressure 

Qf flutter dynamic pressure 

qi t'th generalized coordinate 

a r span of boxes at the rth strip 

V flight velocity 

W work 

W a total work done by system per cycle of oscillation 

Wa work done by strip r per cycle of oscillation 

WJl w r A /w A 

W\ W' A /Sr 

&ij : 3 / 4 slope at 3 /4 -chord point of box i due to deflection of mode j 

u frequency of oscillation 

u; f flutter frequency 
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{*}l/4 


row matrix 
column matrix 
diagonal matrix 
conjugate 
U*] + »[Aj] 

aerodynamic matrix due to rth strip 
imaginary part of aerodynamic matrix 
real part of aerodynamic matrix 

matrix relating boxes to specific strip, defined in equation (15) 
force vector with element i denoting the force acting on box i 
matrix the elements of which are hy at the quarter chord location 
matrix defined by equation (3) 

vector denoting unsteady aerodynamic pressures at the different boxes 
matrix defined by equation (6) 

influence coefficient matrix relating downwash to pressures (see equation (4)) 
vector denoting generalized aerodynamic forces 

vector denoting the generalized complex eigenvector associated with the fluttering mode 
vector denoting the downwash slopes at the 3/4-chord point of the various boxes 
vector denoting deflections at the 1/4-chord point of the various boxes 


GENERAL APPROACH 


The Strip Contribution to the Aerodynamic Matrix 

Assume the lifting surfaces are divided into boxes that form a number of strips along their spans. The doublet lattice 
method (Albano and Rodden, 1969) determines the pressure in each box as a function of the downwash velocities 
{ W / } 3 / 4 at the 3 /4 -chord locations of the different boxes, where 
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where n denotes the number of modes, m denotes the number of boxes, and Oy is the slope of the 3 /4 -chord point 
of box »' as a result of the deformation described by mode j. Similarly, faj describes the deflection at the 3 /4 -chord 
point of box i due to mode j and g,- denotes the ith generalized coordinate. The parameters V and w describe, 
respectively, the fluid velocity and the frequency of oscillation. The deformation can be nondimensionalized 
using a reference semichord length b, that is. 


Therefore, equation (1) assumes the form 



b 


where 


and 


{£} 


r—n 

U> 

II 

an 

021 

ai2 • 
022 • 

■ ain 
• <*2» 

+ ik 

h n 
f 121 

_ ~hn 
f 122 ' 

J • • h\n 

fan 


. &m\ 

Oiml • 

&mn 

3/4 


km2 

h-mn 


, wb 

k= v 

' q i ' 


q 2 

<Zn J 


{?} = 


( 2 ) 


(3) 


4 



The doublet lattice method yields the pressures in each of the boxes in the form 


{p> = [ P]{f } 3/4 (4, 

Substituting equation (3) into equation (4) one obtains 

{p}=lP]{q} (5) 

where 

[P] = [P][Hh/ 4 (6) 

Let the vector {z} j / 4 describe the deflections at the 1 /4 -chord point of each box. The {2)1/4 vector can be written 
in the form 


where 


[ H ] 1/4 = 


Wl/4 

= [Hh /4 {q} 

‘ /*11 

h\2 

■■■ hi n 

f 121 

/ 122 

••• hi „ 

hml 

hm2 

■ ' * h mn 


J 1/4 


(7) 


( 8 ) 


represents the modal deflections at the 1 /4 -chord point of the different boxes. The total force acting on each box is 
therefore given by 


ai 


{F} = 


02 


[ pHq} 


L J 

where o< denotes the area of box i. The virtual work <$( W ) done by these forces is given by 


6(W) = \6q}[H T ] l/4 {F} 


(9) 


so that the generalized forces are given by 


{Q}= [H T ] l/4 {F) 

where [ H t ] denotes the transpose of matrix [H], 

Using equation (9), equation (10) can be written as 


{Q} = [A]{q} 

where 


( 10 ) 


( 11 ) 


[A] = [P r ] 1/4 raJ[P] 


( 12 ) 
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Equation (12) indicates that the element Atj is given by 


m 

= l/4 a k P k j (13) 

*=1 

At this point it is important to note that the summation in equation (13) reflects the contribution to the virtual 
work of all the forces acting on each of the m boxes. Let us consider the contribution to the generalized forces of 
the boxes on the rth strip and denote this contribution by A\j, where 

m 

A r i j = '£h kil/i a k P kJ C kr (14) 

k=l 

where Cjfc r = 0 for all boxes which do not lie on the rth strip and C kr = 1 for all boxes lying on the rth strip. C can 
be arranged in the matrix form 


[C)=[{C k i}{C k2 }---{Cu}} (15) 

where [ C] is of order m x l, with l denoting the number of strips on the lifting surfaces. 

Clearly, A; is obtained by summing up the contributions of all the strips, that is, 

= (10 

r=l 

Based on equations (14) and (15), i aerodynamic matrices [ A] r can be obtained, one for each strip, representing the 
contribution of each strip to the generalized aerodynamic matrix [ A] . Clearly, 

t 

U] = £[A] r (17) 

r=l 


The Work Done by the Aerodynamic Forces 

It is shown in Nissim (1971) that the woric Wa done by the system on its surroundings per one cycle of oscillation 
is given by 


Wa=j L £| (- [Ai + Aj] + i [a r - A t r ] ) (g 0 > (18) 

where {g D } is the complex eigenvector of the fluttering mode and {g*} denotes its complex conjugate. The real 
matrices [Ar] and [ Ai ] are derived from [A] following the definition 

U3 = [A*] + »U/3 (19) 

Following the derivation in the previous section of this work of the strips’ aerodynamic matrices, the work W\ 
done by the rth strip on its surroundings per one cycle of oscillation can be written in a form which is identical to 
equation (18), that is, 
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ff'i = (- [A, + Af]' + . [A R - A$) {,„} 


( 20 ) 



where [ A t Y and [ Ai] r are real matrices defined by 

( 21 ) 


( 22 ) 


= r W\ 

W A =^ (23) 

Sr 

where s r denotes the span of the rth strip. 

It should be noted that the specific energy ratio is the parameter which would indicate the spanwise location 
where most of the energy interchange takes place since it is independent of the strip’s span size. 

Determination of the Spanwise Placement of an Active Control Surface 

The determination of the best spanwise placement of an active control surface for flutter suppression requires the 
following steps to be made: 

Determine the flutter dynamic pressure Qf of the system. 

Increase Qf (so as to lie within the unstable region) and obtain the eigenvector {q 0 } of the unstable mode. 
The amount by which Qp is increased is immaterial since only energy ratios are used herein. In the examples 
treated in this work, Qf is increased so as to lie deep inside the unstable region. 

Compute the matrices [A] and [A] r for all the £ strips, at the reduced frequency appropriate for the unstable 
mode. 

Since the system is unstable, W a should assume a negative value (the system is doing a negative work on its 
surroundings per cycle of oscillation) when using equation (18). 

Use equations (22) and (23) to compute the energy ratios W A and the specific energy ratios W A . Notice that 
negative values of W A and W A indicate that the rth strip is absorbing energy from its surroundings (and thus 
contributing towards instability). 

Determine the location of the strip r = r c where assumes the largest negative value. 

The strip r = r a absorbs most energy per unit span and therefore would present the location where an active 
control surface would be most effective in suppressing flutter. 

The above determination of the most effective placement of an active control surface for flutter suppression, as 
the location where maximum energy per unit span is fed into the system, is a natural consequence of the aerodynamic 
energy approach. It is shown in Nissim (1971; 1977) that a strip that absorbs energy from its surroundings can be 
made, by means of a suitable control law, to essentially stop absorbing energy and even made to be dissipative. On 
this basis, an active control surface placed at the location of maximum energy input will neutralize the maximum 
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3. 

4. 
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[A] T = [A R ] r + i[A I ] T 

At this stage the energy ratio W A of the rth strip will be defined by 


W\ 

W\ = — £- 
A \W A \ 


and the specific energy ratio W A of the rth strip will be defined by 
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energy input and may even turn to dissipate energy. This will clearly lead to maximum effectiveness regarding flutter 
suppression. 

In cases where more than one critical flutter mode exists, the above procedure can be repeated for each flutter 
mode, and the two locations for the control surfaces can be determined. If a compromised single location for the 
control surface cannot be found, then two active control surfaces should be seriously considered. 

It should be noted that 


i 

Wa = Y, W a (24) 

r=l 

Also, since Wa is negative by virtue of equations (2) and (3) above, the use of equations (22) and (24) yields 

= (25) 

r=l 

Note also that the integral over the spans of the specific energy ratios should be equal to — 1 . This follows from 
equation (23), whereby 

W*A = «r w\ 

and following equation (25) 

Ef*3-E*r#'!i-- 1 m 

r=l r= 1 


Example Results 

Three numerical examples are presented with the object of testing the method described above. These examples 
relate to the mathematical model of the drone for aerodynamic and structural testing — advanced research wing 2 
(DAST-ARW2) aircraft used by NASA (Adams and Tiffany, 1985) to test active control concepts for flutter sup- 
pression, and they also relate to the mathematical model of the oblique wing aircraft developed by NASA. Figures 1 , 
2, and 3 show the geometrical descriptions of these models together with the aerodynamic panelling necessary for 
the computation of the aerodynamic coefficients using the doublet lattice method. Figures 1 and 3 also show the 
strip number allocations. In the following, the results obtained are presented and discussed for each of these two 
mathematical models. 

Results for the DAST-ARW2 Mathematical Model 

The DAST-ARW2 mathematical model employs two rigid body modes and ten elastic modes. It yields a flutter 
dynamic pressure Qf = 490 lb/ft 2 (at M = 0 .85) and a flutter frequency up = 117 rad/sec. The unstable eigenvector 
was computed for Q = 550 lb/ft 2 and the matrices [A] and [AY were computed for all 24 strips at the reduced 
frequency k =0.132 associated with the unstable mode. 

Figure 4 shows a plot of the specific energy ratios W A for the 17 strips along the wing. The first four strips 
(not shown in fig. 4) do not exchange any energy since they are on the vertical tail, and the flutter mode involves 
symmetrical mode shapes only. Strips 5 through 21 lie on the wing with strip 5 closest to the root and strip 21 closest 
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to the tip of the wing. It can be seen that the specific energy ratio is negative for all the wing strips, except for 

the root strip where W A is very small and positive. This indicates that all the wing strips (except for the root strip) 
absorb energy from the surroundings, and therefore, all the wing strips contribute to the instability of the system. The 
largest numerical value of W A relates to strip 18, which is the fourth strip inboard of the tip of the wing (coinciding 
with the inboard part of the aileron). Following the method presented herein, this is the best location for the aileron 
for purposes of flutter suppression. For a long aileron (as in the present case), its spanwise placement should be 
such that the aileron midspan point coincides with the midspan of strip 18 (that is, centered around the 80-percent 
span location). The remaining three strips (also not shown in fig. 4) lie on the horizontal tail, and they dissipate a 
negligible amount of eneigy. Note that the integration of W A over the span of the wing and the span of the horizontal 
tail should be equal to — 1 . Note also that the area of W A along the rth strip is equal to the energy ratio W A . In 
conclusion, it can be stated that from the point of view of active flutter suppression, the aileron could be placed a 
little more inboard of its present location. Such an inboard shift can also contribute to an increase in the effectiveness 
of the activated control surface in gust and load alleviation (Nissim and others, 1976). 

Results for the Oblique Wing Mathematical Model 

Figures 2 and 3 show the geometrical description of the oblique wing model with the right wing swept forward at 
an angle of 65° (Burken and others, 1986), together with the doublet lattice aerodynamic panelling and the strip 
number allocations. The model is asymmetric, and it involves five rigid body modes and fifteen elastic modes. A 
flutter computation at Mach 0.95 shows that a mild, 78 rad/sec, vertical tail flutter instability develops around Qf~ 
780 lb/ft 2 . Flutter reaches its peak instability at Q = 1 600 lb/fit 2 and it remains unstable thereafter with the instability 
getting milder as Q increases. This instability can be turned into a neutrally stable oscillation by introducing a 
2.5 percent structural damping. However, even in this latter case, one would expect a large dynamic response of 
this specific mode, with serious effects on the fatigue life of the aircraft, thus warranting the consideration of active 
controls for either stabilization or gust alleviation purposes. The specific energy ratio distribution, as presented in 
this work, was computed for Q = 1600 lb/ft 2 and the results are presented in figure 5. 

It should be noted that strips 1 through 4 lie on the vertical tail, with strip 1 near the root and strip 4 near the 
tip. Strip 5 represents the vertical projection of the fuselage and is not shown in figure 5 (W A in strip 5 equals 
0.011). Strips 6 through 22 lie on the wing, with strip 6 at the left tip of the wing and strip 22 at the right tip of the 
wing. Strips 23 through 26 lie on the left horizontal tail, with strip 23 at the left tip and strip 26 at the root of the 
left horizontal tail. Strips 27 and 28 represent the horizontal projection of the fuselage (from left to right). Finally, 
strips 29 through 32 lie on the right horizontal tail, with strip 29 close to the root and strip 32 at the tip of the right 
horizontal tail. 

Figure 5 shows that most of the energy input into the system takes place through the vertical tail and some 
through the inboard strips of the horizontal tail and the horizontal projection of the fuselage. It should also be noted 
that the wing, the vertical fuselage, and the outboard portions of the horizontal tail play a stabilizing role (with most 
of the wing strips dissipating energy). It is this dissipation of energy that leads to the mild flutter instability exhibited 
in this case. Figure 5 clearly shows that for the suppression of this flutter mode, the active control surface should be 
placed around the center of strip 3 of the vertical tail (that is, around the 60-percent span location). For a long control 
surface, centering the activated control around the 50-percent span location would be more appropriate considering 
the distribution of the specific aerodynamic energy ratios shown in figure 5(a). 

Results for a Modified Oblique Wing Mathematical Model 

The oblique wing model is unique from the point of view of asymmetry. The best placement of an activated control 
surface along its span could present an interesting challenge. Therefore, it was decided to modify the oblique wing 
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mathematical model so as to force the wing to flutter. This was done by considerably lowering one of the torsional 
frequencies of the wing from 45 to 12 Hz. Therefore, the results to be described in the following do not relate to the 
actual wing but to a synthetically modified wing. 

The above modified model of the oblique wing yields two flutter speeds: (1) around Qf = 750 lb/ft 2 and uip = 
79 rad/sec, and (2) around Qp = 1050 lb/ft 2 and up = 70 rad/sec. The first flutter mode is found to be essentially the 
same as the one already described in the previous section of this work (although its mildness is somewhat reduced), 
and, therefore, the results pertaining to this mode will not be repeated. In the following, the results relating to the 
second mode will be described using an eigenvector obtained at Q =1300 lb/fit 2 . 

The specific energy ratio distribution over the different strips of the aircraft is shown in figure 6. It can be seen 
that most of the energy input to the system takes place through the wing strips, with some additional energy input 
through the left horizontal tail and still smaller contributions from the inboard strips of the vertical tail. Some energy 
is dissipated mainly at the outboard portions of the vertical tail, the right horizontal tail, and at some of the wing’s 
strips around its central portion and the extreme right wingtip. Since the wing is essentially responsible for most 
of the energy input to the system, it is interesting to note the effect of the asymmetry of the aircraft on the energy 
input to the wing. Figure 6(b) shows that the strips located around the left tip of the wing are responsible for most 
of the energy input, with the left tip strip (strip 6) showing the largest value for the specific energy ratio. Hence, 
for the suppression of this flutter mode, the activated control surface should be placed as close to the left wingtip as 
is structurally possible. Strips 9 and 10 can also form a reasonable alternative to the extreme tip placement of the 
activated control surface, that is, around the 65-percent span location of the left wing. 

The results presented above regarding the responsibility of the left-most strips of the left wing to the absorption 
of energy which leads to flutter seem surprising and also puzzling. At this stage, it should be stated that the wing 
plane is located 50 in. above the horizontal tail plane. It should be expected that a strong wing-tail interference 
between the left horizontal tail and the tip portion of the left wing will take place. Therefore, it is possible that the 
flutter instability in this case is a result of this interference between the above two lifting surfaces. 

CONCLUDING REMARKS 

The method developed in this work enables the determination of the spanwise placement of activated control 
surfaces for purposes of flutter suppression without resorting to any specific control law. The examples used indicate 
a range of such locations, depending on the flutter mode. It is expected that flutter modes due to external stores will 
indicate still other spanwise placements when compared to those obtained in this work. Where more than one critical 
flutter speed exists requiring neighboring placements of active control surfaces, the placement of a single activated 
control surface at a compromise location can be examined. An example where such a compromise cannot be made 
is given in the modified oblique wing example, where activated controls need to be placed both along the vertical 
tail and the left wing outboard region. 


Ames Research Center 

Dry den Flight Research Facility 

National Aeronautics and Space Administration 

Edwards, California, May 10, 1988 
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Note: 

1. Numbers represent strip locations 

2. VY77X Control surface location 




(b) Horizontal tail planform. 
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Figure 1. DAST-ARW2 model geometric layout, to- 
gether with doublet lattice panelling and strip number 
allocations. 



Figure 2. Oblique wing model — planform view with 
wing skewed 65°. 
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Numbers represent 
strip locations 




(b) Horizontal tail and horizontal projection of the 
fuselage planform. 



(c) Vertical tail sideview. 

Figure 3. Oblique wing model — geometric layout together with doublet lattice panelling and strip numbers. 


13 




When W A >0, strip is 
dissipating energy 
When W A <0, strip is 
absorbing energy 



Figure 4. DAST-ARW2 model, variation of specific energy ratio W a with strip locations along the wing. 
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When W^>0, strip is 
dissipating energy 
When w;<0, strip is 
absorbing energy 
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(b) Wing. 
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(c) Horizontal tail and horizontal fuselage. 


Figure 5. Oblique wing model — variation of specific energy ratio W a with strip locations. 
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When W A >0, strip is 
dissipating energy 
When W A < 0, strip is 
absorbing energy 



(a) Vertical tail. 


When W A >0, strip is 
dissipating energy 
When W A < 0, strip is 
absorbing energy 
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Figure 6. Modified oblique wing model— variation of specific energy ratio W a with strip locations. 
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